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I. INTRODUCTION 

An experiment reported in 1981 [l| measured quan- 
tum noise in a resistively shunted Josephson junction. 
The high frequency noise (both thermal and quantum) 
in the resistor is mixed down to measurable frequencies 
by the non-linearity of the Josephson circuit. This ex- 
periment raised some important questions of principle. 
Quantum noise is not directly measurable, but it niay 
have measurable effects in non-equilibrium situations [2] . 
The experiment has, controversially, been related to the 
question of "dark energy" 0]. The interpretation of the 
experiment involves a quantum Langevin equation (Eq. 
(|1.4p below), and the status of this equation has been 
critically analyzed by one of us Q • The existence of such 
an equation has been demonstrated but only in a simple 
independent-oscillator model 

In the experiment [l[, a resistor R is put in parallel 
with a Josephson junction at temperature T. The criti- 
cal current for the junction is /q and a bias current I > Iq 
is applied. The voltage across the junction is V, and noise 
fluctuations in this voltage are measured. The junction 
has capacitance C. Four different junctions were used 
with slightly different values of the parameters. As an 
example, for junction 2, C = 0.8 pF, Iq ~ 0.5 niA, / = 1 
mA, R = 0.6 11, and one temperature for which details 
are given was T = 4.2 K. The noise is measured at fre- 
quency 183 kHz. 

The theory of the Josephson circuit is particularly sim- 
ple if the capacitance C can be neglected. In [l| a simple 
assumption is made about the dependence on C. One of 
the motivations for this paper is to test the validity of 
this assumption. We find significant deviations for some 
values of the voltage. Our results should be useful if ex- 
periments of a similar kind are performed in the future. 

Our purpose here is to study the solutions of this 
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Langevin equation, and to derive the Josephson fre- 
quency and the generalized impedances and hence predict 
the voltage noise strength . In the interpretation of the 
experiment [l[, certain approximations were made, and 
we examine the validity of these approximations. We 
use two forms of perturbation theory, and also numerical 
methods. 

A mechanical model for the dynamics of the Langevin 
equation (equation ()1.4|) below) has been studied exper- 
imentally and theoretically in [6|. 

The fundamental theoretical input Q is the quantum 
version of the Nyquist fluctuation-dissipation theorem 
(originally derived in Q). In the present context, it re- 
lates the current noise in the resistor at temperature T 
to the conductance 1/R, and states 



1 



where 



/2ft^l^ 




'hn' 




coth 






2T 



[Siin) + Sii-n)] = — coth — (1.1) 



Siin)5in - n') = (/„(r!)/„(-r!')> (i-2) 

/„(r2) being the Fourier transform of the current noise 
/„(i) in the resistor, and the expectation value (...) is 
with respect to a thermal distribution at temperature T 
(where we use units such that Boltzmann's constant is 

The right hand side of (|l.ip has the finite limit 
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(1.3) 
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as T — 0, and this offers the possibility to measure quan- 
tum noise. For this it is essential to have the symmetrized 
version of Sj; the unsymmetrized Sj{^l) itself has a van- 
ishing low-temperature limit for f2 > 0. This reflects the 
impossibility of getting energy directly out of the vac- 
uum. It is, of course, essential that /„ is a quantum, 
non-commuting, operator. 

The right hand side of (|l.ip contains the distribution 
function for a quantum oscillator of frequency fl at tem- 
perature T. However, as emphasized in 0, this by no 



2 



means implies that the resistor contains oscillators, or 
even that the system is bosonic. 

The second theoretical input is the assumed quantum 
Langevin equation, connecting the voltage and voltage 
noise to the current and current noise in the shunted 
Josephson circuit: 
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In sin 9 = I + In 



(1.4) 



where 9 is the Josephson phase (including noise) and C 
is the capacitance of the junction. Since /„ is a quantum 
operator, 9 must be also. The status of this equation is 
less secure than that of (see Q)- The voltage V 
across the junction is given by 



Vit) 



h 



(1.5) 



We will always assume that I > I^. Then, in the absence 
of the noise /„, the solutions 9c{t — t) of (|1.4p settle down, 
after the decay of transients, to a periodic dependence on 
time, with frequency ilj = 27r/r,/, so that 



9,{t + Tj - t) = 9c{t - t) + 27r 



(1.6) 



for all t. The arbitrary constant t is included in order to 
make explicit the infinite set of solution depending on the 
initial condition (although (jl.4p is a second order equa- 
tion, the initial value of 6 is not an independent initial 
condition for the steady periodic solution). 

On the right of (|1.4p , the current noise /„ is certainly a 
quantum operator, but / is presumably a classical quan- 
tity. Since /„ is small, we can write 

9 = 9c + 9n (1.7) 

where 9c is the (classical) solution and 6'„ is the small 
quantum noise. Then the latter obeys, to a good ap- 
proximation, the linear quantum equation 
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in which cos 9c is a c-numbcr coefficient. 

Because of (|1.6p , there is a Fourier series (we use j for 
^/^) 

cos9c{t-i) = ^gkexp{jknjt)exp{-jkmj) 

k 

= Y,ak{i)exp{jknjt) (1.9) 

k 

with k an integer and g'^ = g^k- We define a hybrid 
Fourier integral/series by 

Init)=Y, dnink{^)cx^[jt{n + k^lj)] (1.10) 



and similarly for 9 and V = {h/2e)d9/dt. Then the 
transform of (jl.Sp is 

I]fc'[-P'fe(^)4,fe' + h9k-k'{i)]0nk' 

= Ek' Xk,k'm0nk'm = Ink{^)- (1.11) 



where 



Fkin) = - (^^^ {n + kQjf 



Eq. ()1.11|) has solution 

9nkm^j2^kk'minkm 



(1.13) 



where the matrix Y is the inverse of X (for each il). 
We assume that there is a non-zero limit 



so that (by (|1.5 



^im [jf]yo,fe(f^)] = ^^fe. 



Vno{0)=J2ZkInkiO). 



(1.14) 



(1.15) 



This is certainly true for C = [T3|, and also consistent 
with our conclusions in the remainder of this paper. 

We have not made explicit the t dependence oi X, Y 
and Z, but it follows from (HH]), ([TTU)) . ([LT^ and (|ni)) 
that Z is proportional to cxp{~ j ktQ j) . 

We define a periodic function /„ (t) associated with /„ 
by 



init) = ^/„fe(r! = o)exp{jtknj) 



(1.16) 



and similarly for Vn Then, for the purpose of computing 
the Zk defined in (|1.14l) it is sufficient to use /„ and Vn- 

The basic problem is to obtain from p. II) and (|1.8p 
information about the voltage noise (using (jl.Sp ). One 
approach to this problem 0, [Tlj is to simulate the current 
noise numerically, but it is not clear how to do this consis- 
tently with the frequency dependence in p.ip . Also, the 
noise is supposed to be a quantum operator. Our method 
is indirect. We first obtain the impedances Z^ defined in 
(|1.15p . Since these are assumed to be classical quanti- 
ties (this is the consequence of the assumed existence of 
a quantum Langevin equation with classical coefficients 
in it), to find them it is sufficient to take the noise in 
(|1.8p to be classical. Having found the impedances Zk, 
we can use (jl.l5p (with quantum noise) together with 
(|l.ip to determine the strength of the (quantum) volt- 
age noise. An advantage of this method is that (|l.ip is 
exactly respected. 

Thus, the purpose of this paper is to study the solu- 
tions of equation (|1.4p and (jl.8p , and the calculation of 
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rij and Zk- This enables us to check the vahdity of some 
of the approximations made in [l|. In section Hill we use 
perturbation theory for small C. This turns out to be a 
series in C^, and we are able to calculate the O(C^) con- 
tribution to f2j, and to Zk but for fc > 1 only. In section 
\V\ we use perturbation theory in l/I. Again, there are 
only even terms, but wc calculate 0(1//^) contributions 
to Qj and to Zk- But this approximation is not useful for 
the values of / in the experiments. In section lVTl we give 
results of numerical calculations of flj, Ro and Zi,Z2, 
for several values of /3 and i (dimensionless parameters 
defined in (|2.ip and (|2jp below). We differ from some 
other numerical work j7l. in that we do not attempt 
to simulate the actual noise (which should be consistent 
with (jl.ip ). but just determine the Z^ factors. Our con- 
clusions are given in graphs, which may be useful in the 
interpretation of experiments. 

One particular conclusion is that, in the range < /3 < 
0.2 (with /3 defined as in (|2.ip below) Zi decreases and 
Rd increases. In Table |ll] we show by what factors the 
values of Z^ and of {Zi/RdY are predicted to change 
between /3 = and /3 = 0.38. As we shall discuss, these 
ratios are new corrections to the predicted voltage noise. 

We conclude that the apparent quantitative agreement 
between theory and observation in [l| may not be as good 
as it appears. Certainly, our calculations should be rele- 
vant to any future repetition of this type of experiment. 
For another possible applications of equation <\\A\ . see 
Rcf. [H. 



where in general ujj is the natural frequency of the solu- 
tions of ()2.2|) (when transient decaying terms have died 
out), and ujjq is the value in this approximation (/3 ~ 0): 



(2.4) 



Similarly, denotes the solution of (|2.2p when = 0. In 
this case, equation (jl.l4p gives only three non-zero values 
of Zfc, with k = 0,±1. 

The dynamic resistance is defined to be 



Rd = -rr, td 
clI 



dv 

ZQ = -77, 

at 



(2.5) 



where V denotes the time-average over one cycle. In this 
approximation (/3 = 0) 



27r 

TJO 



(2.6) 



so, from ((^ . 



(2.7) 



II. NOTATION AND A SIMPLE SOLUTION 



Equation (jl.isp gives 



It is useful to define the dimensionless variables 
/ V 2eRIo 



2eR'CIo h Zk 
, uj = \l, Zk — — . 2.1 



In terms of these variables, and denoting dd/dr by 0', 
equation (jl.4p becomes in the absence of the current noise 



pe" + e' + : 



(2.2) 



In the special case when j3 is negligible, this equation 
is easily soluble. A particularly convenient form of a 
solution is 
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sm Wo = 

cos Of) = 

dBQ_ ^ 
dT 



1 -t- i sin(ti;,7o''') 
i + sin(a;joT) 
LP JO cos{ujjot) 
i + sin(a;,/oT) ' 

i + sin(a;,7oT) ' 



(2.3) 



{Vnomonmt - 5]cxpHif],,(A:-fc')] 



fc,fc' 



X ZkZl.il^kiOKk-m (2.8) 



where we have made explicit the i dependence in (|1.14p 
inherited from (|1.9p . We now average (|2.8p over t through 
one period Tj and use p.2p to get 



SviO) = (Ko(0)y;o(0)) 



di{Vno{Q)v:,mi 



1 



(2.9) 



Since Z^ = Z^k, the right hand side of this equation 
automatically contains the symmetrized products of cur- 
rents, as in equation (jl.ip . Finally using (fTTTI) with 
= k^lj, (|2.9p gives the required prediction for the low 
frequency voltage noise S'y(O). 
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k\zk\^ujj coth(pfcwj 



(2.10) 



k>l 



where we have used p.ip and p = eRIo/T (and p is close 
to 1 in the experiment. 

For the special case C ~ 0, wc have 



RI 



Z±i\ = 



Rh 



and all other Zk^O. Then (|2TT0)) gives 
ATRP eR^ll 



(2.11) 
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/2- 
X coth 
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(2.12) 



Using (P?T]) . (1131), (^31) and (P^ . this may aUcrnatively 
be written 
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coth 



(2.13) 



It has been suggested [J that in this form, if observed 
values of V and Rd are used, it may also be a good 
approximation for non-zero /3 (defined in (|2.ip ). One 
purpose of this paper is to test this approximation. We 
note from (l2l0l) that the quantity Sv{'^)/R\, in (l2J3l) 
depends upon the magnitudes of the ratios {Zk/Ro)'^- 



III. PERTURBATION THEORY FOR SMALL 
CAPACITANCE 

This section is about the calculation of the Joscphson 
frequency; it is not concerned with noise, so all variables 
are classical ones. 

Since /3 in equation (|2.2p has modest values in the ex- 
periments [l| , one might expect perturbation theory in /3 
to be useful. This is especially so as the series turns out 
to be an expansion in /J^. This property follows because 
is invariant under 



T->-T, 13^-13, e^n^e, i^+i. (3.1) 
To construct the perturbation series, we expand 

= 0n + 1301 +(3^02 + ..., (3.2) 

and 

TJ = TJO + pTjl + ^V,72 + . . . , (3.3) 

where r = 27r/cjj is the Josephson period and tjq is given 
by (1131). Inserting into 



0[ + (cos6lo)6li 



^0 ' 



(3.4) 



0^ + (cos0o)02 = -0i+ysin0o, (3.5) 

where prime denotes differentiation with respect to r. 
Equation (|3.4p gives, using 



[i + sin(wjoT)6'i]' = - 



wj^Q cos(a;jor) 



Hence 



i + sin(a;joT) 
i + sin(w,7o''') 



(3.6) 



(3.7) 



where we have chosen the initial value 0{Q) = This 
solution has the same frequency wjo as 6*0; so ri = in 
p.3p . This is in accordance with our remark above that 
the expansion is series in /S^. 

From the definition of T, for any r, to the requisite 
order, 

2tT = 0{t + Tj) - e{T) 

= 0O{.T + Tj) ~ 0o{t) 

+ P^[e2{T + Tjo)~6{T% (3.8) 



givmg 



TJ2 = -[02(t + T,o) - 02{T)]le'j^{T). (3.9) 



(To be consistent, the right hand side of this equation 
must turn out to be independent of r.) 

Then, from dlSl), and dSH), and using (g^l). 



TJ2 



T + Tjo 



dr'li + sin(cL'jor')] 



X [0'l{T')-l{sin0o)0,{r'f] 
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r+Tjo 



dr' 



sin(wjor') In |^1 + 

1 + i sin(a;jor') 

(i + sin(a;jo'r'))^ 

sin(a;joT') 



sin(w,7oT') 



X InM 1 



(3.10) 



where we have integrated by parts the 0'{ term in the first 
line. It turns out that all the integrals arising in p.lOp 
which contain a power of ln(l-|-sin((u;joT)/z) can be done 
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by integration by parts: 



TJ2 



Where 
Hir) 



/T+Tja 
dr' 

wjo cos^(aj,/oT') dH 
{i + sm{LOjoT'))^ dr' 

cos{ujjot) \ 1 



(3.11) 



i + sin(a;joT) 



■ InM 1 



smujjoT 



In 1 



sin(cL'joT 



(3.12) 



It is now a simple matter to eomplete the integration and 
find 



TJ2 = -^jo 



r+Tja 



dr 



{i + sin(wjor'))2 
= -2tt{i~ujjo). (3.13) 

Thus the Josephson period to second order is 

Tj = Tjo[l - fiojoit - wjo)], (3.14) 

and the frequency is 

= ^,7o[(l + P'^^joii - ^,7o)]- (3.15) 



From (j2.6p . in this approximation the dynamical resis- 
tance is 



ro = zo = V P [i- w,7o)(2i - wjo). (3.16) 

^,70 



Figures [2] and [3] show a comparison of p.lSp and p.l6p 
with numerical results. 



IV. PERTURBATION THEORY FOR 
IMPEDANCE OF HIGHER HARMONICS 

In this section, we use perturbation theory in j3 to 
study the impedances z^. Because zi has a contribution 
(see (|2.1ip ) to zeroth order in /3, 

|zi|^ = |zio + /3zii+/32zi2 + ...p 

and we need the second order term because of its in- 
terference with the zeroth order one; so first order per- 
turbation theory on its own has no physical significance. 
We have found second order perturbation theory pro- 
hibitively complicated. But for > 1, there is no zeroth 
order contribution and so first order perturbation theory 
is relevant by itself, and this we now study. 

Since our aim is to calculate the impedances Zk (which 
are classical quantities), it is sufficient to introduce, in- 
stead of the quantum noise 7„, a classical driving force at 



the Josephson frequency wj into (|2.2p . To our approxi- 
mation, it is sufficient to use Eq. (2.4). Thus we use the 
equation 

pe" + e' + sine = i + ak exp{jkujjoT) (4.1) 
where ak is complex and infinitesimal. We expand 

9 = 00+ ak9a, + bOi + bauOia^ + ... (4.2) 

Then 

^'ak + ^ifc cos 6*0 ^ cxp{ikujjo) (4.3) 

so 

[{i + sm{ujjoT))]ea, = F 
= d-^'li^ + sin(wjoT')] cxpijkujjoT'). (4.4) 



(4.5) 



6i is given in (|3.7p . 
To the next order. 

So, using and 
d 



U! 



dr 

2 1 -I- i sin(wjoT) 
i + sin(ajjotT) 



[(i + sin(w,/ot))6'iaJ ^ 

sm{ujjoT) 



In 



- [i + sin(wjoT)] 6'^^. 

Wc want to determine the period of 9 in 
period be 



Tj = TJO + bakTla,. + 



then 



27r = 9o{t + T JO + bakTla^) - doi^r) 

+ bak[9iakiT + Tjq) - ^la^^{T)] 

where r is arbitrary, so 

riaAir) = -[0iaAr + rjo)-9,aAr)]- 
It follows from (lL4l) and (l46l) that 



(4.6) 
Let this 

(4.7) 



(4.8) 



(4.9) 



T + Tjo 



dr' 



2 1 -f i sin(a;jor') 



<J0 



{i + sin(a;joT'))2 



1 n , smwjoT 
X ln[l + r ]9a,iT) 

i^ + siniu;Jor'))9';Jr')] 



(4.10) 



Integrating by parts and using the periodicity of the in- 
tegrand, the second term in the square bracket in (|4.10p 
is converted to 



UJjQ Sm{uJjoT')9ia^. 



(4.11) 
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Then, using (|4.10|) gives 



wjoJr I sin(a;joT') 



1 + ln(l 



sm[u)jQT 



-G'(r') (4.12) 



where 



i + sin(a;,7o''') 
1 + In I 1 



sin(a;,70'r) 



(4.13) 



Again, because of the periodicity, the contribution from 
the G' term in (|4.12p is zero. Finally, inserting the value 
of F' from (|4.4[) and omitting an integral of a differential 
of a periodic function, we get 

Y PT + Tjo 

Tiofc = / dr' cos(wjot') 

WJO Jt 



X In 



1 



sin(w7oT') 



exp(jfcwjor'). (4.14) 



Because the integrand in (|4.14[) is periodic, we may 
replace the limits of integration by and tjq Integrating 
by parts again, it gives 



1 r™^,/ cosjujjoT') 
2wjo Jo i + sin(wjoT') 



exp{j{k + l)ujjotT') cxp{j(k - l)ojjnT') 



1 



1 



(4.15) 



The integral may now be done by the substitution C = 
exp(ja; jqt') and integrating round the unit circle. There 
are poles at C = —ju, —ju~^ where 



u = 1 + 0; JO, u 



(4.16) 



and —ju ^ lies within the unit circle (taking i > 1 for 
the moment). Remembering that 



the result is 



-fe+i 



-fe-i 



k - 1 



1 



(4.17) 



(4.18) 



(The restriction i ^ I may now be relaxed.) 

For low enough temperature, or high enough current /, 
when p in (2.10) is considerably greater than 1, equation 
(|2.10P for the voltage noise approximately contains the 
sum 



fc>i 



Zk\ 



(4.19) 



which, in the small (3 approximation (|4.18|) , may be eval- 
uated to give 
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i + LO.J0 



(4.20) 



In the experiment [l|, fc is of order 1, but still (|4.19p is 
a lower bound to the contributions from fc > 1 (since the 
coth factors are greater than 1). 



V. PERTURBATION THEORY FOR HIGH 
CURRENT 

In this section, we develop an expansion in powers of 
1/i = /o //. This is algebraically simpler than the expan- 
sion in /3, but is relevant only to a small range of / in the 
experiments. 

Because of the invariance of (|2.2p under i —i,9 
—6, the expansion is in powers of l/i^. 

Define s = ir and a function i?(s) — 0{t). Then (|4.ip 
becomes 



+ + - sin ^ 
i 



1 + — exp(jfcs), (5.1) 



and we will be concerned only with fc = 1 or 2 in this 
section. In general the frequency of the last term in (|5.ip 
ought to be a; J, but to our order it is sufficient to use the 
zcroth approximation wjo ~ i. Since the driving force is 
classical, is also classical in this section and in section 

EI. 

Expand 



^ = 'do + -^i + ^^2 
, "1 q , "1 q 



la2 



(5.2) 



with similar expansions of tj and wj. (We will find that 
Tnak = for 71 < fc + 1.) We take t?o = s, being the solu- 
tion without a transient decreasing exponential. Then 



(i/3)?9'/ + = - sin?9o = - sin s 



(5.3) 



(z/3)<,^+^;„^ =exp(js), (5.4) 
(i/3)?9^' + 19^ = -(cosz?o)i?i = -(coss)t9i, (5.5) 

(«/3)^L + ^2a, = -(cOSZ?o)l^lai = - (cOS s)l?i,, . (5.6) 



The solution of (|5.3p with no transient, and a conve- 
nient initial value, is 



exp(js) exp(-js) 



where 



1 + j5 1 - 29 



9 = H 



(5.7) 



(5.8) 
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Since (|5.7p is periodic with the same period as t^o, the 
period is not changed to first order, as expected. 
The solution of (lOl) is 



-jexp(js) 



(5.9) 



Again, the period is unahered, and so there is no correc- 
tion to V or to zq at this order. 
Inserting (|5.7I) into 



(i/3)i9 



1 



4(1+52) 
X [2 + (l-j5)exp(2js) 
+ (l+jg)exp(-2j.s)]. 



so 



(5.10) 



(5.11) 



It foUows, in a similar manner to (|4.8p . the period at this 
order is 



1 



2tt 



and 



TJO + ^TJ2 = — 



UJj = I 



1 



1 



2^2(1 + ^2) 



(5.12) 



1 - 



dojj 
di 



1 

2^2(1 + g2) 
2^2(1+^2)2 



(5.13) 



to this order. (For /3 = 0, this is consistent with (|2.4p ). 
Inserting (|5.9p into (|5.6p . we find that 



hence 



exp(2js) 



[2{j ~ g){2j - 4g)] [2(j-.g)]' 
1 

■&2a2{s + 2TT)~da2{s) = -TT- 



Therefore 

and 

zi 



T2ai = TTyr 



(5.14) 



(5.15) 



(5.16) 
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;T2ai 



2i^ 



(5.17) 

Comparisons of (|5.13|) and ()5.17p with numerical data 
are shown in Figs. [H [3] and SI The perturbation theory 
is in reasonable agreement for i > 2.0. 

We can use the expansion in 1/i for k = 2, but we 
must go to one order higher in 1 /i to find Z2 ■ We use the 
02 terms in (|5.ip and (|5.2p . The equations are 



ffC+^U = exp(2js), 

5^?2a.+^2a. = -(cOSs)^i„,, (5.18) 
5^^302+^302 = -(cOSs)7?2a2 + (Sins)7?it9ia2. 



Then 



^2a, = 



exp(2js) 
2 j - 4g 

exp(js) 



[2(i-5)(2i-4g)] 
exp(3js) 
[2(95 -3j)(4g- 2.7)]^ 



(5.19) 



where the second term will turn out to be irrelevant. 
Putting (f5TT9| into the last of ([5T8| and using ([577] 



5^^302 +^^302 



1 



[4(j-g)(2j-45)] 
1 

[4(2i-45)(j+5)] 



(5.20) 



where the omitted terms are periodic and irrelevant. 
Hence 



Z?3a2(s + 27r)-t93a2(s) 



27r5 



[4(j- 25)(52 + l)]' 
Thus we find that, to this order, 

\9\ 



\Z2\ 



[4*2(52 + 1)^(452 + 1)] 



(5.21) 



(5.22) 



As expected, this is zero \i (3 — Q (since 5 = i/3). 

The above formulae (|5.13p . (j5.17p and (j5.22p are com- 
pared with numerical results in Figs. [2l [3l |4] and [6l re- 
spectively. As expected, the perturbation theory is not 
accurate for the smaller values of i, but it agrees quite 
well with the numerical results for i > 2. 



VI. NUMERICAL RESULTS 



We start with equation 
terms: 



including real forcing 



pe" + 0' + sin6' = i + as sin(a;'T) + cos{ijj't), (6.1) 

where Og and Oc are small (in practise we have verified 
that, in the range 10~^ < Os, Oc < 10~^, linearity holds 
within the numerical precision of the calculations), and 
uj' is adjusted so that 



wj(iu', as, ac) = w'. 



(6.2) 



LOj being the frequency of the solutions of (|6.ip (after the 
decay of transients). 

The numerical solution of Eq. (jO.ip was investigated 
using several initial conditions for both i < 1 and i > 1. 
Our main interest here are the cases when i > 1. In these 
cases, we start with arbitrary initial conditions at r = 0, 
and begin measurements on the solution at r = tq where 
To is about 100 (when transients have become negligibly 
small, for all relevant choices of /3, i and tu'). Then, after 
N = 10^ cycles, the condition 6'(ti) - 9{tq) = 27rA^ yields 
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1.09705 



(Oj 1.09695 



1.09685 



1.09685 



1.09695 

to' 



1.09705 




0.2 0.4 0.6 0.1 
P 



1.2 1.4 



FIG. 1: Graphical representation of the equation 
ujj{ijj',as,0) = ij' ■ The crosses are the numerical re- 
sults for ujj{uj',as,0) where as = 10~^, 2 x 10"'', ... and 
7 X 10"". In this example /3 = 0.45 and i = 1.4142. 




0.2 0.4 0.6 0.! 

P 



1.4 



FIG. 2: The dimensionless average voltage as a function of /3 
and four values of i. The continuous and broken lines are the 
perturbative results in powers of 1/i and /? respectively. 



FIG. 3: The dimensionless dynamical resistance as a function 
of P and four values of i. The continuous and broken lines are 
the perturbative results in powers of 1/i and /3 respectively. 



p 


4^ X 100 

V 


^ X 100 

rn 


0.00 


6.066 


11.611 


0.30 


6.923 


13.380 


0.60 


4.371 


12.790 


0.90 


1.921 


7.760 


1.20 


0.854 


4.168 


1.50 


0.415 


2.264 


1.80 


0.221 


1.289 


2.10 


0.126 


0.772 


2.40 


0.077 


0.485 


2.70 


0.049 


0.318 


3.00 


0.033 


0.216 


3.30 


0.023 


0.152 


3.60 


0.016 


0.109 



TABLE I: Differences between the numerical and the 1/i per- 
turbation theory for both and ro, for i = 1.41. 



the period tj = (ri — to)/N (the size of the time step 
employed in the numerical code was At = 10~^, between 
T = and T = tq, and At = 10~® after the transients). 
At the solutions of (|6.2p . but not otherwise, the results 
for ojj = 2tt/tj and \z\ are stable against changes in the 
initial conditions and in tq ~ 100 and TV ~ lO'^. 

The next step in the numerical procedure consists in 
solving (|6.2|) . In the figure [1] we illustrate this procedure 
with an example of the intersections of a^, 0) with 

the diagonal line for /3 = 0.45, i = 1.4142 and a series of 
values of as . Although figure [1] only shows a few points 
in the vicinity of the diagonal line, the numerical code 
employed 40 times more points, in the same range shown 
in figure [1] in order to reach the required precision in 
the determination of the impedances. In this range of 
variation of Os , the point of intersection increases linearly 
with tts- 

In general the numerical procedure solves (j6.2p and 



yields the result 

iAjj{as,ac) = wj(0, 0) + asZs + a^Zc + ■ ■ ■ (6.3) 

which gives the Josephson frequency and the two real 
impedances (for k ^ 1). (We also have obtained the 
nonlinear terms in (|6.3p which are not important in the 
present analysis). Using this approach, we were able to 
perform a detailed numerical calculation of the quantities 
wj, ro, \zi\ = i/zf + zf and \zi\/ro- We also have in- 
vestigated the higher harmonics (A: = 2, 3, ... ) and com- 
puted the corresponding quantities such as |z2|, \z2\lr1:) 
and |zip-)-2|z2p. Our main interest is to investigate how 
significant is the dependence on /? for several values of i. 

Let us start with the results for v = uj[Q,Q) and ro 
(see Eq. (|5.13p ). These are shown in figures [2] and [3] as 
functions of p and four values of i. From these figures 
one can see what are the numerical values of /? and i such 
that the perturbative results, obtained in the previous 
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FIG. 4: The impedance associated with the first harmonic as 
a function of /3 and four values of i. The full lines represent 
the 1/i perturbative result. 



FIG. 7: \zif + 21^2 P as a function of /? and four values of i. 
The full lines are the 1/i perturbative results. 
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FIG. 5: Same as in figure|4]but divided by ro- 



FIG. 8: Same as in figure [7] but divided by td- 
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FIG. 6: The impedance associated with the second harmonic 
as a function of /? and four values of i. The full lines are the 
1/i perturbative results and the dashed lines are the order P 
perturbative result given by (|4.18|l . 



sections, can be trusted. For instance, the perturbative 



increase with 13 only persists up to about /3 = 0.3, for 
i = 1.4; after that rjy decreases. Although the numerical 
computation has been performed up to /3 = 4, here we 
are focusing on the interval which may be more realistic 
from the phenomenological point of view. Besides, larger 
values of (3 are very well described by the 1/i perturbation 
theory even for i = 1.4. For instance, we have found 
that, for i = 1.4, the difference between the numerical 
and the 1/i perturbative result for r^) is smaller than 1% 
for (3 > 2.0 (see table IJ). 

Let us now consider the results for the impedances as- 
sociated with the forcing terms. The numerical results 
for \zi\ ~ y/zf+z^ and \zi\/rD are shown in figures H] 
andO respectively. Similarly to the case of the dynam- 
ical resistance the curves for the 1/i perturbative result 
(full lines) underestimate the exact numerical points by 
an amount which becomes negligible when i increases. 
Also, the perturbative result becomes indistinguishable 
from the exact result for /? sufficiently large. 

The results for the impedances associated with second 
harmonic are shown in figure [5] The validity of the per- 
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i 


ri 


r2 


1.41 


0.96 


0.90 


1.91 


0.75 


0.72 


2.41 


0.55 


0.56 


2.91 


0.44 


0.49 
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TABLE II: The ratios n = [Zi{l3 = 0.38)/Zi(/J = 0)]^ and 
ra = liZi/RD)(P = 0.38)/(Zi/i?B)(/3 = 0)]^ are shown in 
columns 2 and 3, respectively, for the values of i, shown in 
column 1. Notice that these ratios are the same if we use 
either Zi and Rd or zi and ro- 



turbative results for small /3 or large i can be seem in 
the figure. As in the previous cases the curves for the 
l/i perturbativc result become very close to the exact 
numerical result for f3 sufficiently large. Finally, in fig- 
ures [7] and [5] we shown the results for |zip + 2|z2p and 
(I Zip + 2|z2p)/r£). These quantities are relevant for the 
calculation of the low frequency voltage noise (according 
to the Eqs. (jTTU)) and glUl)). 



VII. CONCLUSIONS 



We have obtained the generalized impedances Zk for 
a range of values of the two dimensionless parameters 13 
and i. These allow one to find the current noise strength 
from (2.10). We are particularly concerned with the de- 
pendence on the parameter (3 which measures the impor- 
tance of the capacitance of the junction. For /? = the 
Langevin equation (1.4) has a simple analytic solution. 
An important question is how far the results for /? = 
provide a good approximation when /? is of order 0.5. 
Figures 2 to 8 and table [Til exhibit the dependence on /3, 
which leads to significant corrections to the voltage noise. 

Since the analysis used in [l[ neglected the /3 depen- 
dence [l^ the conclusions of that paper may not be as 
decisive as claimed. If similar experiments are done in 
the future, our results should be useful. 
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